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Abstract 

We present a method for constructing superoscillatory functions the superoscillatory part of 
which approximates a given polynomial with arbitrarily small error in a fixed interval. These 
functions are obtained as the product of the polynomial with a sufficiently flat, bandlimited envelope 
function whose Fourier transform has at least N — 1 continuous derivatives and an A^-th derivative 
of bounded variation, N being the order of the polynomial. Polynomials of arbitrarily high order 
can be approximated if the Fourier transform of the envelope is smooth, i.e. a bump function. 
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I. INTRODUCTION 


Superoscillation is the counter-intuitive property of a bandlimited function to oscillate 
with local frequencies that are larger than its maximum frequency component. Early ex¬ 
plicit reports of this property are found in the held of signals and systems BEl. although 
controlling the oscillations of bandlimited functions through their zeros is a much older sub¬ 
ject m- From the viewpoint of physics, superoscillations imply the ability of a radiating or 
imaging system to produce or resolve wave features much hner than the bandwidth of the 
system suggests. Elements of the phenomenon can therefore be traced back in the classical 
quest for superdirective antennas (originating from Einstein’s concept of needle radiation) 
[1] and for optical imaging beyond the diffraction-limit [^. Interestingly, the term super¬ 
oscillation was actually born much later within quantum mechanics as part of the concept of 
weak measurements EE!, which subsequently motivated the systematic analysis of faster- 
than-Fourier functions m and their temporal or spatial evolution as quantum or optical 
wavefunctions 110]. In recent years, superoscillations enjoyed a renewed interest within the 
held of super-resolution imaging, as optical technology enabled their use for focusing light in 
the far held at spatial scales below the dihraction limit mm- In any application involving 
superoscillations, however, one has to hght against the inevitable diminishing of the signal 
amplitude in the superoscillatory region, as follows from standard Fourier analysis na. 

Physical implications aside, the construction of superoscillatory functions is an interest¬ 
ing mathematical problem in itself. Two general approaches can be distinguished in the 
literature so far. The hrst uses the sampling theorem to express a superoscillatory function 
as a series of shifted sine functions with only a hnite number of non-zero coefficients. The 
latter are determined through a linear system of amplitude constraints on a hnite grid of 
points mm- The grid is hner than a Nyquist sampling grid in order to make the function 
super-oscillate. This approach is also related to the old problem of constructing bandlimited 
functions with a given set of zeros |3]. For this particular problem, there is also the possi¬ 
bility to directly replace some of the zeros of a bandlimited function with the desired ones 
without changing its bandwidth. The zero-replacement theorem has appeared in several 
independent works [21 ESI, as well as in explicit connection with superoscillations HU. 

The second approach involves expressing the superoscillatory function as a Fourier inte¬ 
gral and again imposing a set of amplitude constraints on a hne grid of points. Variational 
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techniques are subsequently used to find the minimum-energy function that solves the prob¬ 
lem [131 HE]- Constraints for the derivative of the function can be applied as well [13]. The 
variational approach too has old roots in the field of information theory [2D]. 

In addition to the above methods, the literature provides well-studied examples of super- 
oscillatory functions, such as the familiar f(t) = (cost -|- iasint)^ 0 HD], as well as specific 
Fourier-type integral representations where superoscillations emerge due to the presence in 
the transform of a delta-function-like factor centered in the complex plane [S]. 

A common feature of the existing general methods is that the superoscillatory function 
is required to satisfy a discrete set of constraints, concerning the value of the function or its 
derivative on a fine grid of points. This however implies limited control over the actual shape 
of the function in the superoscillatory interval. In sampling methods in particular, the points 
at which the sine basis functions are centered can be chosen quite arbitrarily na leading 
to an infinite number of different functions satisfying the same superoscillatory constraints. 
A similar lack of control over the shape of superoscillations exists when specific methods or 
prototype functions are used. Their parameters can be tuned to define the maximum local 
frequency or the total number of the superoscillations but there is little or no flexibility at 
all in controlling their shape. 

An immediate question is whether one can construct superoscillatory functions in a more 
continuous way, namely having control over the shape of the function at least in the super¬ 
oscillatory region. More specifically, one may ask if it is possible to construct a superoscil¬ 
latory function f{t) G L‘^{R) (square integrable in the entire real line) that approximates a 
desired analytic function p{t) with arbitrary accuracy in a certain interval {—a, a). The ac¬ 
curacy of the approximation can be quantified in terms of a norm of the difference f{t) —p{t) 
over the interval, say ||/ — ( 7112 , namely 

/ ra \ 1/2 

yj [f{t) -P{t)f dtj <e, (1) 

for some small e. Of course, for an arbitrary analytic function p{t), one can generally speak 
only of an approximation because a strict equality f(t) = p(t) over a finite interval would 
imply, by analyticity, an equality for all t which is generally impossible since p{t) may not 
even belong in L‘^{R), as for example in the case of a polynomial or a trigonometric function. 

The answer to the above question is positive. To prove this, we here present a simple 
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method that allows to construct superoscillatory functions that approximate a given polyno¬ 
mial PN{t) {N being the order) with arbitrarily small error in a hnite interval. Although the 
method refers to polynomials, general analytic functions p{t) can be treated too by first ex¬ 
panding them into their Taylor series around some to ^ ®)- The expansion is truncated 

to some order N so that the corresponding Taylor polynomial pN{t) is an approximation of 
pit) in this interval in the sense \\p — Pn \\2 < ci- Then a superoscillatory function /(t) can 
be found that approximates this polynomial in the sense ||/ — PAr ||2 < ^ 2 - By the triangle 
inequality ||/ —p ||2 < \ \f ~Pn \\2 +\ \p — Pn \\2 < ^1 + ^2 which is a number that can be made 
smaller than any given e. 


II. METHOD 

We consider the Paley-Wiener space PBA of real square integrable functions /(t) G 
whose Fourier transform P(co') is supported on [—7r,7r], namely bandlimited functions with 
hnite energy and bandwidth tt. We also assume the real polynomial 

N 

PN{t) = ^ant'^ ( 2 ) 

n=0 

as the target or desired shape of our function in the superoscillatory interval {—a, a). For 
having at least two zeros in this interval, one requires N >2. Now consider a known function 
e(t) G PWt, that we term the envelope function. We additionally assume that its Fourier 
transform Eiu) is at least C^~^{—oo, 00 ) (i.e. it has at least N — 1 continuous derivatives 
for all real u) and has a N-th derivative of bounded variation. Then by the familiar property 
of Fourier analysis the function fit) = PAr(t)e(t) has the Fourier transform 

N 

= (3) 

n=0 

where the superscript (n) indicates the n-th derivative with respect to the argument. By 
our assumptions, the above transform is clearly zero for |a;| > tt and of bounded variation 
(hence integrable, p > 0), therefore fit) G PBA-. 

The above can be stated alternatively using the smoothness-decay property of Fourier 
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transforms [21]. If the Fourier transform of the envelope function is at least oo, oo) 

and has an A^-th derivative of bounded variation, then e{t) decays at least as as 

|t| —>■ oo, so that p]\f(t)e(t) is still a function in PW.„, despite the polynomial growth of 
PN{t). 

Moreover, if e(0) = 1 and e(t) is sufficiently flat at t = 0, one has f(t) ^ PN(t) in the 
interval of interest {—a, a). The flatness of e{t) can be independently controlled through a 
dilation transformation e{t) —)■ e{t/D) where D > 1. Our superoscillatory function hnally 
reads 


f{t) = pN{t) e (4) 

Note that the dilation conhnes the spectrum of e{t) even more hence e{t/D) is still a member 
of PW^. 


III. AN EXAMPLE 


As an example consider the function 


f{t) =P3{t) e 



3^3 /P 

~ir 



(5) 


where sinc(a;) = sin(7ra;)/(7ra;). It is easy to show, by successive convolutions of the Fourier 
transform of the sine function (which is equal to 1 for |a;| < tt and zero otherwise), that the 
function e{t) = sinc"^(t/4) has the Fourier transform 


E(i^) 


16 

/ e{t)e-^^^dt = — 

J —oo ^ 


a;|^ - TTu'^ + 

I (tt - |a;|)^ 


u 



I < |a;| < TT 


( 6 ) 


and zero for |a;| > vr (plotted in Fig. [^c)). The above is C'^(—oo, oo) differentiable with an 
integrable third derivative hence, according to the previous discussion, sinc^(t/4) is an ap¬ 
propriate envelope to multiply with a cubic polynomial and obtain a /(t) G PW-,,. For s < 1 
and H > 1, the f{t) of Eq. ([^ superoscillates in (—s, s) approximately as the cubic polyno- 
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FIG. 1. (a) The function of Eq. (§ for s = 0.1 and D = 1. The inset zooms into the superoscillatory 
interval (—0.1,0.!), also showing the corresponding polynomial P 3 (t) in red. (b) The imaginary 
part of the Fourier transform of f{t) (the real part is zero), (c) The Fourier transform of the 
envelope function (d) The error ||/ — P 3 II 2 in (—0.1,0.!) as a function of the dilation factor D. 


mial Psit) = “ V®]- The prefactor sets the amplitude of the superoscillation to 

1 . 


An example of the function of Eq. ([^ is shown in Fig. [^a) for s = 0.1 and D = 1. 
Notice that its Fourier transform, shown in Fig. Bb). is hardly distinguishable from the 
third derivative of the Fourier transform of sinc^(f/4). Indeed, according to Eq. (|^, we 
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have for our example 


which shows that the third derivative term dominates for s << 1. However, even in the 
case of vanishingly small s, the hrst derivative term cannot be neglected as it corresponds 
to the linear term in the target polynomial Psit) which is responsible for its oscillatory (or 
sinusoidal) shape. This verihes the nature of superoscillations as a delicate interference effect, 
something that has already been noted in the literature [10] from a different perspective (our 
present point of view is that of the spectrum). 


IV. ON THE SPECTRUM 

Generalizing the above remark, consider a function f{t) G PW.,r that superoscillates in 
(—a, a) with the oscillation amplitude in the order of unity. In the context of our approach, 
the behaviour of f{t) in (—a, a) can be approximated with the help of a polynomial as for ex¬ 
ample its truncated Taylor series around t = 0. If the superoscillations have a characteristic 
scale s << 1, the Taylor polynomial will be of the form 

N 

PN(t) = 

n=0 

with the real coefficients being in the order of 1. Stating that alternatively, the magnitude 
of the derivative /^'^^(O) is in the order of s“"’. This is easily understood by considering 
polynomials with roots separated by multiples of s, for example 

In the context of our method, the entire superoscillatory function f{t) is the product of such 
an approximating polynomial with an envelope function e{t) that decays sufficiently fast as 
|t| —)■ cxD and is flat enough around t = 0 with e(0) = 1. Multiplying Eq. (|^ with such an 
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envelope and taking the Fourier transform we obtain 


^ -n 
n=0 ^ 

Thus F{u) contains contributions of the derivatives of E{uj) weighted with s~'^ which makes 
the contribution of the highest-order derivative dominant. Nevertheless, despite the 

spectrum of /(f) resembling its superoscillatory behaviour is actually due to the 

progressively weaker, essentially perturbative, contributions of the lower-order derivatives 
(cn) with n < N. This explains why it is difficult to tell if a bandlimited function 
is superoscillatory from its spectrum (“there is no hint of superoscillations in the power 
spectrum” HU]). 

The smallness of the amplitude of superoscillations can also be quantihed in the context 
of our approach. Having normalized their amplitude to unity, the amplitude of the function 
increases outside the interval (—a, a) according to the approximating polynomial pwit). The 
maximum is reached at some t ^ D and, according to Eq. (§, is in the order of (D/s)^. 
Since N determines the number of zeros in the superoscillatory interval, we have verihed the 
exponential dependence of the amplitude on the number of superoscillations |13j . 


V. CONCLUSION 

We have presented a simple method for designing superoscillatory functions that ap¬ 
proximate a given polynomial PAr(f) with arbitrarily low error within a given interval. The 
functions are obtained by multiplying the polynomial with a sufficiently flat bandlimited 
envelope function e{t) whose Fourier transform is (at least) as smooth as needed to keep 
the product in the Payley-Wiener space. Appropriate envelope functions can be obtained, 
for example, as powers of the sine function or, more generally, as the Fourier transform of 
polynomial splines. Envelope functions that work with polynomials of arbitrarily high order 

N can also be constructed. Such are the inverse Fourier transforms of bump functions, 

— 1 

namely C'°°(—oo, oo) functions with a compact support, as for example E{uj) = e for 

|a;| < 71 and E{u) = 0 for |a;| > vr. We have hnally seen that the superoscillatory nature 
of a bandlimited function constructed with present method is hidden in its spectrum as a 
series of geometrically diminishing contributions from derivatives of the envelope’s Fourier 



transform. 
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